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I. INTRODUCTION 



Parity-violating (PV) inelastic electron-nucleus scattering is an important tool in the 
study of hadron structure [Q. In combination with PV elastic electron-proton (e-p) scatter- 
ing, measurements of the PV quasielastic (QE) electron-deuteron (e-d) asymmetry allow a 
separate determination of the strangeness magnetic form factor, G« (Q 2 ), and the isovector 
axial vector form factor, Cr^ T (Q 2 ). Knowledge of Gm\Q 2 ) provides a window on the 
role played by sea-quarks in the electromagnetic structure of the nucleon. The axial vector 
form factor, in contrast, is sensitive to nucleon structure effects in higher-order, electroweak 
radiative corrections. These corrections, which depend on the species of lepton probe (hence, 
the "e" superscript), share features with corrections relevant to other precision electroweak 
measurements, such as the PV asymmetry in polarized neutron /3-decay. The proper inter- 
pretation of such measurements relies on an adequate understanding of electroweak radiative 
corrections [§]. 

Recently, the SAMPLE collaboration has performed a separate determination of G$ 
and Gi e){T=1) at Q 2 = 0.1 (GeV/c) 2 using PV e-p and PV QE e-d scattering [|| |, |. 
The results indicate a value for G$ < ' T_1 ' ) consistent with zero. At tree-level, one expects 
Ga T (Q 2 = 0) = —1.267, while radiative corrections reduce the magnitude by roughly 
40 ±20% H [Tj. These corrections include potentially significant hadronic contributions that 
are responsible for the estimated theoretical uncertainty. To make the measured value of 
Gi e)(T=1) close to zero would require additional effects not included in the calculation of 
Refs. § 0. 

One possibility, which we explore in this paper, is the contribution from the PV nucleon- 
nucleon (NN) interaction. The latter induces small admixtures of opposite parity states into 
the deuteron as well as the scattered np partial waves. These parity admixtures contribute 
to the PV asymmetry when a 7 is exchanged between the electron and target 1 . Moreover, 
in contrast to the effect of Z° exchange, these hadronic PV effects in 7-exchange give rise to 
a term in the asymmetry which does not vanish at Q 2 = 0. For sufficiently small Q 2 , this 
term would dominate the asymmetry. One might ask, then, whether the omission of this 
term in the interpretation of the SAMPLE deuterium measurement is responsible for the 

1 Current conservation also implies the presence of PV electromagnetic two-body current contributions, as 
we discuss below. 
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apparent, anomalously large radiative corrections entering G^ . 

Below, we show that the magnitude of this Q 2 -independent hadronic PV contribution is 
too small to account for the observed ^ ^ effect. Based on simple scaling arguments, the 
relative importance of the Q 2 -independent contribution - compared to the "canonical" Z°- 
exchange induced asymmetry - goes as ~ 10~ 4 m^/Q 2 . Thus, at the SAMPLE kinematics, 
Q 2 = 0.1 (GeV/c) 2 , we expect the hadronic PV contribution to generate at most a few parts 
in a thousand correction to the asymmetry - far short of what would be needed to close the 
gap between the theoretical and experimental values for ^ T . 

We also carry out an explicit calculation of the hadronic PV contribution and verify 
the expectations based on these scaling arguments. Our computation follows on the work 
of Refs. [§, H, in which the hadronic PV contribution to PV threshold deuteron electro- 
disintegration was studied, and the calculation of Ref. [[R]], which treated PV QE e-d 
scattering. In the latter analysis, only parity-mixing in the deuteron wave function was 
considered. In the present study, we also include the contributions from parity mixing 
in the final e-d scattering states as well as from PV two-body currents. Our results are 
consistent with both of these earlier calculations, but give a more complete treatment of the 
QE case. 

The remainder of the paper is organized as follows. In Section ||, we review the for- 
malism for PV QE scattering and hadronic PV, identify the relevant operators and matrix 
elements to be computed, and present the scaling arguments for the relative magnitude for 
the hadronic PV contribution. Section |T| gives a discussion of the bound and scattering 
state wave functions, which we determine first in the plane wave approximation and subse- 
quently using the Argonne Vig potential. We present the results of our calculation in Section 
pV| , where we consider two cases: threshold electro-disintegration and QE scattering. Figs. 
|7|, H, and [L3|, which show various contributions to the PV asymmetries as a function of 
Q 2 , summarize the main results of this work. A summary discussion appears in Section [V]. 
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II. PV ELECTRON SCATTERING AND HADRONIC PV 



A. Basic Formalism 



The PV asymmetry for inclusive e-d scattering of an unpolarized target can be expressed 
in terms of two response functions: W EM , the parity-conserving (PC) electromagnetic (EM) 
response, and W pv , the PV response arising from the interference of EM and PV neu- 
tral current amplitudes. One may decompose the former in terms of the longitudinal and 
transverse response functions 



wEM = EW^+^w] > w 

w=Ef—Ei 
f 1 

fKq) = J2 F cM, (2) 

J>0 

F 2 T {q) = J2i F EM + KA<l)} , (3) 
j>i 

where vl,t are the standard kinematic coefficients (will be defined later), q^ = (u,q) is 
the four momentum transfered into the nuclear system (Ei and Ef are its initial and final 
energy). The F X j(q), X = C,E, M, are the charge, transverse electric, and transverse mag- 
netic multipole matrix elements depending on q — \q\. They are defined through multipole 
operators, 6 C = M,6 E = f e \ and 6 M = if ma 9 0, [11], p|, as 2 

Fxj ^ = ~7W^f Yl (-1) T/_Mt ( T/ T ^ \ (Jf, T /-'-Oj T (qy:':Ji,T^ , (4) 
V2^ + 1 T t^i I -M T; M Ti / 



where the :: denotes reduced matrix elements in both angular momentum and isospin 0. In 
the spherical basis, while a collective quantum label a refers to (E a , L a , S a , J a , Mj a ,T a , M Ta ), 
the sum runs over all indexes except Ef and Mj f because they have been carried out to 
get Eq. ®. 



2 The extra "i" for T ma9 is introduced so that when real wave functions are used, the form factor Fm is 
real. 
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For the PV response, one has 

W PV(Z) 



E 

/ 



v L Wj( v (q) + v T Wl v {q) + v T iW^' A (q) 
W^ v {q) = -g e A J2 F cAl)Fcj(q), 

J>0 



LU = Ef 



J>1 



W? A (q) = -g e v J2 [ F ej(q)FmjM +F M j(q)F Ej5 (q) 



(5) 
(6) 
(7) 
(8) 



where the Fx (5) refer to weak, neutral current multipole matrix elements and the "5" subscript 
indicates multipole projections of the axial vector current. The Fx (5) are defined in a similar 
fashion as Eq. @ - up to different coupling constants, however, for the axial form factors, 
it is M 5 and T eU ° which have additional factors of % while T ma95 is without one [||] 3 . The 
kinematic coefficients, vl, vt, and vt> are 

(Q 2 /q 2 )\ (9) 
(Q 2 /g 2 ) 2 /2 + tan 2 (£/2), (10) 
y / (g 2 /g 2 ) 2 + tan 2 (0/2) tan(0/2) , (11) 

where Q 2 = q 2 — u 2 , 9 is the scattering angle of electron. 

The PV response functions W A y arise from electron axial vector (A) x hadronic vector 
current (V) interactions, while Wy A is generated by the V(e) x A(had.) interaction. At 
tree-level in the Standard Model (SM), the electron vector coupling to the Z° is suppressed: 
g v = — 1 + 4 sin 2 9w ~ — 0.1 (the axial vector coupling is g e A = 1). 

In terms of these response functions, the PV QE asymmetry due to Z°- exchanges is 



vl 
vt 



A 



(z) 



GaQ 2 W PV(Z) 



For quasielastic kinematics, u and q are related, viz, u ps q 2 /2m N . 



(12) 



B. Hadronic PV Effects 

Hadronic PV effects in the target generate 0(a) corrections to the tree-level contributions 
for Wy' A . These hadronic PV effects arise when a photon, rather than a Z°, is exchanged 

3 In Ref. Fx j and Fxj 5 are defined with an extra 1/2 and —1/2 factor. It is found that the minus sign 
in Ref. is a typographical error. The 1/2 factor is absorbed here in the overall scale A® R . 
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Figure 1: Contributions due to hadronic PV in elastic e-N scattering, where (g> denotes the PV 
meson-nucleon coupling. 

between the electron and hadron. Because the vector -fee coupling Q e = — 1 is an order of 
magnitude larger than g*, one expects the relative importance of the hadronic PV effects - 
compared to the tree-level amplitude - to be of order 

8V2ira 1 



9tt 



-0.01 



(13) 



G^A 2 , 1 - 4 sin 2 B w g A 

where g n = 3.8 x 10~ 8 sets the scale for hadronic PV interactions, A x = 4irF n « 1 GeV 
gives the scale of chiral symmetry breaking, and g A = 1.267 ± 0.004. In the case of elastic 
e-N scattering, hadronic PV arises via diagrams of the type in Fig. [I]. These corrections 
induce a PV jNN coupling, or anapole moment. The latter has been used in the one-body 
estimate of R A given in Refs. [§, ^j. Those analyses indicate hadronic PV-induced anapole 
corrections of —6 ± 20%. 

Two-body hadronic PV contributions arise from the diagrams in Fig. ||. Figs. ||(a) 
and (b) indicate parity-mixing in the initial and final state wave functions, while Fig. ||(c) 
indicates the PV two-body EM current contribution. Each contributes to an effective axial 
vector EM transition amplitude, whose effects appear as corrections to the F X j s (q) multipole 
matrix elements. 

In computing the parity-mixing matrix elements, we use the model PV Hamiltonian given 
in Ref. @: 
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Figure 2: Contributions due to two-body hadronic PV in e-d scattering. Here, M and M' denote 
the identities of mesons. 
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H PV (r) = 9 * N ^ hlT (in x t 2 )z{<ti + ^2) ■ ^(r) - 
4v2m JV 



x [(1 + //„)zct! x a 2 ■ u p (r) + {B x - a 2 ) ■ v p (r)] 
9.nn ( h o + ^ + fa) 



277^ \ u 2 

x [(1 + /j s )z<7! x <r 2 ] ■ Uu{r) + (cti - B 2 ) ■ v u (r)] + 

- -7—(n ~ T 2 ) z (ai + B 2 ) ■ {g^NNhlv^ir) - g P NNh l p v p {f)) 

- ^^'(^1 >< ^1 + **) ■ u p {f) , (14) 

(x) 

where the h M and c/mnn are the weak, PV and strong, PC meson-nucleon couplings, 

respectively; and u a (f) = [pi — p 2 , e~ mar /r], v a {r) = {p± — p 2 , e~ mo,r /r}. Values for the 

(x) 

h M have been predicted theoretically using a variety of approaches. For purposes of our 
calculation, we will adopt the so-called DDH "best values" and "reasonable ranges" of Ref. 



13fl . The latter are consistent with constraints obtained from a variety of hadronic and 



nuclear PV experiments, as discussed in Refs. [[y|, [T^ . 

Current conservation requires that one include contributions from PV two-body currents 
to the PV multipole matrix elements. These currents have been derived in Ref. |16| , [i~7|| from 
the diagrams in Fig. ||(c). Complete expressions for the coordinate space current operators 
could be found in Ref. [[[7]]. To illustrate the structure of these operators, however, we 
give the complete two-body current operator associated with the ^-exchange component 
of H PV : 

J^x u x 2 Y v - = ~ e9 * NNh * {T X ■ r 2 - t 1z t 2z ) (^(y- xt) + a 2 5®(y - x 2 ) 
8v27rm Ar \ 



-k*i ■ Vi - B 2 ■ V 2 )(f + \^-V y ){5^\y- x x ) + 6®(y- x 2 ))) 



x 

\2 



/i = 1,2,3 



« 0(v/c) 2 /2 = } (15) 

where Xi denotes the position of the 7-th nucleon and x = \x\ = \x x — x 2 \. Inclusion of 
these PV two-body currents guarantees that the the continuity equation, V • J = i[H,p], is 
satisfied at the operator level of the PV NN interaction. 
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Current conservation may also be implemented in writing down the various multipole 
operators. Since nuclear model calculations based on realistic potentials generally break 
current conservation, it is useful to implement the latter via the form of the multipole 
operator. A well-known example is Siegert's theorem |I8], which allows one to rewrite the 
J = 1 transverse electric multipole operator, Tf =l , in terms of the electric dipole operator 
(the J = 1 charge multiple) in the long wavelength limit. An extended version of this 
theorem [|l9|1 allows one to implement the constraints of current conservation for transverse 



electric operators of arbitrary J and momentum transfer. In general, one has [20 



T?(q) = S e J(q) + E;J( q ). (16) 



For J = 1, one has 



Sf(q) = ^j- J d 3 x[p(x)x,H], (17) 

Rfii) = ~jj d " x iffiW • f x lift > ( 18 ) 

where is the vector spherical harmonic. 

Note that in the long wavelength limit, Sf{q) gives the leading contribution to Feu Feu 
and Fei 5 - For elastic electron scattering, hermiticity and time-reversal invariance require 
that Fei and F E \ must vanish. Moreover, contributions from Sf 5 (q) also vanish for elastic 
scattering, since the commutator in Eq. fllTD leads to a factor of uj = Ef — Ei = 0. A non- 
vanishing contribution arises from _Rf 5 (g), whose matrix element constitutes the nuclear 
anapole moment contribution. The latter is proportional to Q 2 which cancels the l/Q 2 from 
the photon propagator to produce a Q 2 -independent scattering amplitude at lowest order. 
This contribution is kinematically indistinguishable from the Z°-exchange contribution to 
Feu an d; thus, represents a simple multiplicative correction to the tree-level axial vector 
response. 

For the inelastic transition of interest here, matrix elements of Sf 5 {q) do not vanish, nor 
do they contain a power of Q 2 to cancel the l/Q 2 from the photon propagator (q 2 w Q 2 
at low energies). The scattering amplitude associated with this operator goes as l/Q 2 at 
low-Q 2 . When multiplied by the factor of Q 2 in Eq. (|T2j), this term thus generates a Q 2 - 
independent, non- vanishing contribution to A LR at low-Q 2 , in contrast to the Z°-exchange 
asymmetry which vanishes at Q 2 = 0. For sufficiently low-Q 2 , the Q 2 -independent hadronic 
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PV EM contribution will dominate the asymmetry. As we show below, such a kinematic 
region may be reached in principle with threshold PV electro-disintegration. For PV QE 
scattering at the kinematics of Ref. ||, however, the effects of hadronic PV appear to be 
negligible. 

At the kinematics of QE electron scattering, a tower of final state partial waves contribute 
to the amplitude, and one must include the effects of J > 1 multipole matrix elements 
(both transverse electric and transverse magnetic). This situation contrasts with threshold 
disintegration, where only the lowest J partial waves may be reached. As we show in Section 



rVj , the multipole contributions to the PV QE asymmetry (due to hadronic PV) saturate 
for J ~ 7. All multipole matrix elements having J > 1 carry factors of Q 2 , so that they do 
not contribute to the Q 2 -independent term in the PV asymmetry. Nevertheless, the sum of 
their effects represents a tiny correction to the Z°-exchange asymmetry. 

It is useful to illustrate how PV NN effects contribute to the various multipole matrix 
elements entering the axial response, Wy' A . Consider, for example, a transition from the 
deuteron ground state to the 1 S , continuum state. In the absence of the tensor force com- 
ponent of the strong NN interaction, the deuteron ground state is pure 3 Si. The PV NN 
interaction will mix P-states into these S'-waves. In ordinary perturbation theory, one has 



| 3 5x) -> l 3 ^) + |35 a ) , (19) 
I 1 S ) -> \ 1 S ) + \ { S ), (20) 



where the parity mixtures (denoted with a tilde " ") are 



to = e i v'T'i'^ - < 2i » 

fc=l,3 

I'So) = \ 3 P°} {3P i":i?° } - (22) 

For this Jj = 1 to Jf — transition, only J = 1 multipole operators contribute. For the 
PC 7-exchange contribution, one has only the magnetic dipole transition between the un- 
mixed 3 5*i and 1 S initial and final state components, resulting in a non-zero Fmi form factor. 
For the PV Z°-exchange amplitude, only the operator Tf 5 (q) contributes, connecting the 
un-mixed 3 Si and 1 Sq components and leading to a non-zero Fei 5 - The PV NN interaction 
also contributes to the latter in three ways: (a) a non- vanishing matrix element of Tf(q) 
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between the initial state \ 3 S\) and the final state \ 3 Po) parity admixture, (b) non-vanishing 
matrix elements of Tf(q) between the | 1,3 Pi) mixture in the initial state and the final state 
l 1 ^); (c) matrix elements of the PV two-body current operator Tf 5 (q) between the | 3 £i) 
and \ 1 Sq) components. All other contributions are higher-order in the weak interaction and 
can be neglected. The analysis is similar when the D-state components of the deuteron and 
scattering state induced by the tensor force are included, as is the analysis for transitions to 
higher partial waves. 

Because the hadronic PV contributes to the asymmetry by inducing axial photonic cou- 
plings, to incorporate these in the expression of Eq. ([12]), one only has to modify axial form 
factors Fxj 5 by 

F X j 5 -> I' X.I:, + (31-x.h, ■ 

where 



F XJs ^F X j B +j3m, (23) 



P G nQ 2 9% ' { } 



The EM axial form factors may be decomposed as 



T f T Ti 



= y v ., 

V2 Ti + 1 ^ x \ -M Tf M Ti 

u f /r f \\()^iq)\\j^r,) + u' f :f f \\()J :i ( q )\\j,.r,) + u f /if\()^(q)\\J,r,) k?5) 



X 

where |Jj,Tj), (Jf,Tf\, and Oj^(q) (the two-body PV EM operators) represent the effects 
caused by hadronic PV. One may then express the asymmetry due to hadronic PV (through 
the radiative corrections) as 

A^ R = 2 4r 7 , (26) 



where 



(27) 

u)=Et-Ei 



f j>i 

A simple scaling argument allows us to estimate the relative impact of the two-body 
hadronic PV contribution. For backward-angle scattering as studied in the SAMPLE exper- 
iments, vt> ~ vt ^> vl, so the ratio of asymmetry due to hadronic PV, Eq. (|26|), and Z°, 
Eq. (|12|) is 

4« 8V2na W PV ™ 8V2na (j], v ) 



(28) 
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While (j z ) at backward angles is dominated by the magnetic NC component and scales as 
(a), (jp V ) scales as (a) times an additional factor introduced by hadronic PV. Using Eq. 
(Pf) for a guidance, this factor is roughly 

VA 7 ' 1 (29) 



8y/2ir \ m N x 

With g n NN — 13.45, and h n ~ 4.5 x 1CT 7 , the scaling rule is 

4(7) ^2 /„-m n x\ 

where we have also included the NC magnetic form factor Gm{Q 2 ) in the denominator. For 
small Q 2 , one has Gm ~ Hv = 4.70. Taking x ~ 1 fm, (e~ m7TX / (m N x)) ~ 0.1, therefore at 
g 2 ~ 0.1 (GeV/c) 2 , we have A[ 7 i/4fl in the order of a few 0.1%. 



III. TWO-BODY WAVE FUNCTIONS 

In order to compute PV matrix elements in Eq. (f25|), we need two-body wave functions. 
The latter are solutions of the Schrodinger equation 

(ff + ff PV )|V + $ = £|V + $. (31) 
Since H PV is much smaller than H , first-order perturbation should work well in this process. 
That is, Eq.(|3l"D can be solved in two steps: first, the PC part \ip) is determined by solving 

(E-H )\4) = 0, (32) 

and second, the PV admixture \ip) is determined from 

{E - H )$) = H PV \i/,) , (33) 

with \ip) obtained in the first step. In what follows, we explore two different approaches, 
one using the plane wave approximation - which ignores the final state strong interaction 
(FSI) - and one using a potential model calculation, which includes the FSI. 



A. Plane Wave Approximation 

Although the plane wave approximation is naive and simple, we employ it as a toy-model 



calculation to achieve some initial insights. In addition, the computation of Ref. (10| em 



ployed a plane wave Green's function to compute the PV admixture in the deuteron, though 
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un-mixed deuteron wavefunction was obtained using the Bonn potential. By comparing the 
plane wave computation with the potential model solution (see below), we hope to obtain a 
sense of the errors introduced by the plane wave approximation. 

In this approach, all the radial components of scattering partial waves are spherical Bessel 
function, j'i(pr), where L is the relative orbital angular momentum and p is the relative 
momentum. 

The parity admixture, in first-order perturbation expansion, is expressed as 

l^> = E^)p J -p-<^W, (34) 

where \<f)) forms a complete eigenbasis. This could be computed if one knows the Green's 
function 

G (^=£<M<f>. (35) 

In the plane wave basis, closed-form Green's functions exist for both calculations of deuteron 
and final state mixing. 

For the deuteron mixing, the transition involves a bound state (binding energy Eb < 0) 
to continuum state transition, therefore, E B — E^ < 0. The Green's function is 

G^(x,y) = J2 

L,S,J,Mj,T,M T \ 71 J 

x ^( 7 r<)^( 7 r>)^S(^)^(^)xi MT X^, (36) 



where 7 = ^ym N \E B \; iL(l r ) and k L (^r) are the modified spherical Bessel functions of the 
first and third kind; r< (r>) refers to the smaller (larger) radial coordinate of x and y; y 
and x denote the spin-angular and isospin wave functions. The factor Slst, which is 1 if 
L + S + T is an odd number and otherwise, enforces the generalized Pauli principle. 

For the final state mixing, because of the pole at = E^, we must add a small imaginary 
number ±ie to the energy denominator, as in the scattering problem. In this way, we obtain a 
retarded (advanced) Green's function corresponding to the — ie (-He) prescription. However, 
only the real part of this Green's function gives a non- vanishing response function. The real 
part is equivalent to the average of retarded and advanced ones 

G ( ^(x,y) = hsr(m N p T ) 

L,S,J,Mj,T,Mt 

x jL(^r<)n L (^r>)^(^)^(^)xr T XT T , (37) 
13 



where p F is the relative momentum of the final state and ni(p :F r) is the spherical Neumann 
wave functions. 

Illustrative results for the plane wave calculation are given in Figs. [3| and §. We note 
that, in comparison with the complete, coupled channel potential model computation (see 
below), use of the plane wave Green's function overestimates the degree of parity-mixing 
in the deuteron ground state. For parity-mixing in the scattering states, we also find a 
mismatch between the two approaches, though no systematic pattern emerges as to the 
magnitude or sign of the difference. The problem may be particularly severe for the 3 Si and 
3 Di scattering states which, in the plane wave approach, are not automatically orthogonal 
to the deuteron wave function \T>). Although one might attempt to solve this problem by 
implement orthogonality by hand, viz, 



it is questionable whether this ad hoc solution is rigorously correct. For these reasons, then, 
we rely only on the coupled channel potential model computation to determine the nuclear 
PV contribution to the inelastic asymmetry. 

B. Potential Model Calculation 

Although the NN potential determined directly from solving QCD is not available, a 
variety of modern phenomenological potentials successfully fit NN scattering data (below 
350 MeV or so) as well as deuteron properties with reasonable x 2 values. Here, we use the 
Argonne Vis potential (AVis) pl]1 - 

The PC scattering wave function, 



where u{r) denotes the radial wave function; k = ^/m N E; and the overall constant is fixed 
by the normalization condition (E'...\E...) = 5(E' — E)..., is determined by solving the 
Schrodinger equation. This task is eventually reduced to integrating a one-dimensional 
differential equation for the radial component and solving for the phase shift. However, due 
to the tensor force, for J > 0, states having quantum numbers (L, S, J) = (J — 1, 1, J) and 
(J + 1, 1, J) are coupled, requiring that one solve a coupled set of differential equations. 




3 S 1 , 3 D 1 )-\V)(V\ 3 S 1 , 3 D 1 ) 



(38) 




(39) 
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The normalization for radial wave functions are fixed by their asymptotic forms. For the 
uncoupled channel problem, one has 



u JLS (r -»• oo) = r sin(fcr + Ltt/2 + 5 JLS ) , (40) 

where 5 JLS denotes the phase shift. For the coupled channel problem, the convention intro- 
duced by Blatt and Biedenharn (BB) [22|], with two eigenphases shifts 5^}\ b^f 1 and a mixing 
parameter e J? is adopted. 4 The two orthogonal, real solutions are 

uilj : \ / cose 7 sm(kr + (J- 1)tt/2 + <5? } ) \ 

m = r J 1 1 ^ ' J , (41) 

42 J+1 / \ v sin ej sin(A;r + (J+l)7r/2 + 5 ( / 1) ) J 

wSj-! V , /-sine J sin(fcr + (J-l)7r/2 + 5 ( / 2) )\ 

(r -> oo) = r . (42) 

42./ +1 / \ cose, sin(A;r + (J + l)7r/2 + 5 ( ; ) ) y 

It should be noted that while we will still refer to solution 1(2) as 3 [J — 1] j( 3 [J + 1] j) state, 
it contains a component involving the other channel. We have verified our calculations by 
reproducing the experimental phase shifts. 
The deuteron wave function, 



(f\v, Mj) = I ~Y- y& + y% | x" , (43) 

is obtained by solving the eigenvalue problem for binding energy Eb and D/S ratio. The 
asymptotic and normalization conditions are 

u(r <C 1) oc r , u(r ^> 1) oc r k (jr) (44) 

u>(r <C 1) oc r 3 , u>(r ^> 1) oc r k 2 (jr) (45) 

dr [w 2 (r) + w 2 (r)] = 1. (46) 



Although one can follow a similar strategy and obtain the PV wave functions by the 
Green's function method mentioned in previous subsection, it is not straightforward to do 
so; the unperturbed wave functions are too complex to allow one to obtain analytical results 



4 The "nuclear bar" convention, defined in Ref. |2^|, is more commonly used in literature. However, the 
phase parameters of these two conventions are totally interchangeable. The reason for our choice is that 
all the wave functions are purely real in BB convention. 
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as in the case of plane waves. Therefore, following the same approach as in Ref. we 
directly solve the inhomogeneous equation, Eq. fl33|) . 

The basic idea is to solve the problem twice, once with the "source" term off (thus a 
homogeneous equation as solving the PC wave function) and then with the source term 
on. A general solution for the inhomogeneous equation, ip g , can be expressed as a linear 
combination of solutions for the homogeneous equation, called the complimentary solutions, 
ip c (i), plus the particular solution, ip. Therefore, in order to obtain the particular solution, 



the complimentary part has to be fully subtracted. Thus, we must determine a«, i = 1...N, 
N being the number of coupled equations. 

In the case of solving scattering wave function, the asymptotic behaviors of both ip g and 
ip c can be expressed as linear combinations of incoming and outgoing spherical waves. While 
the interactions cause phase shifts of outgoing waves, the incoming waves are not altered. 
This observation tells us that, ip, the parity-mixed component induced by the PV NN 
interaction, should not contain any incoming component. Using this result, a{ are solution 
when the incoming wave components of ip g and ^ c (i) completely cancel in Eq.([47|). 

Except for 1 Sq and 3 Pq, which can only be mixed into each other, all the other uncoupled 
states, 1,3 L J=L could have mixtures from 3 [L — 1} L and 3 [L + 1} L states. For the coupled 
states, 3 [L = J — 1] j and 3 [L = J + l]j, both mix to 1,3 Jj. If the mixture is an uncoupled 
state, and we have 




(47) 



a e 



i(kr—Lir) < t i(kr—Ln) 



(48) 
(49) 



ce 



i(kr—L-rr) < j i(kr-Lir) 



then 



(50) 



When the mixtures are coupled, a two channel calculation is needed. If one has 



^(i = l,2 




o(z, l) e - i ( fcr -( J - 1 W + b(i, i) e i ( fcr -( J - 1 )^) 
a(i, 2)e" i ( fcr -( J+1 W + 6(i, 2)e i( - kr ^ J+1 ^ 

C(l, l) e -i(kr-(J~l)Tr) + d ^ ^ e i(kr-(J-l)n) 

c(i, 2) e - i ( fcr -( J+1 W + d(i, 2)e i ( /cr -( J+1 ^ 




(51) 



(52) 
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then 



/ oti(i) i 









^ = 4(1) - ai(l)4(l) - a 2 (l)^(2) , (53) 
= 4(2) - «i(2)4(l) - a 2 (2)&(2) , (54) 

with 

'a(l,l) a(2,l)\ 1 / c(z,l)\ 
a(l,2) a(2,2) J ^ c(i, 2) J ' 
Note that all the coefficients here are complex. This implies that the mixed wave functions 
are also complex. However, in our framework, only the real part will contribute to the 
response function W /py(7) and, thus, to the asymmetry. 

Various criteria exist for testing the numerical solutions: i) they should satisfy the differ- 
ential equation, ii) they should be independent of the initial conditions used to integrate the 
differential equation, iii) they should be proportional to the source term, i.e., if the source 
term doubles, the solution should also double. These conditions are employed to make sure 
we obtain the correct solutions. 

As for the parity admixture of deuteron, it is determined using the same procedure. 
Since one is dealing with a bound state however, the asymptotic behavior is given by a 
linear combination of modified spherical Bessel functions, ii and k^. The physically realistic 
solution is obtained by completely subtracting the ii component, because it exponentially 
diverges. 



IV. RESULTS AND DISCUSSION 



First, we compare the two approaches discussed in Section [TTJ. Fig. |3] shows an example 



of the comparison between the plane wave scattering states and the ones obtained from 
the potential model calculations (in our case, it is AVis). Though at E rd = 1 MeV, the 
3 Po solutions look almost the same, the plane wave 1 S , state differs from the more realistic 
solution by a large phase shift as well as in its radial shape at small distances. Note that the 
latter difference is important because the PV NN interaction is very sensitive to the short 
range behavior of wave functions. Therefore, the plane wave approximation is not adequate. 
Fig. H shows the PV mixtures for these 1 S and 3 P states, 3 P an d 1 S respectively. They 
are similar to the results of Ref. [§], which were obtained by using Reid soft-core potential, 
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10 20 30 40 50 

r(fm) 



10 20 30 40 50 

r(fm) 



(a) 1 5o channel (b) 3 Pi channel 

Figure 3: Comparison of scattering state wave functions: dashed lines give the plane wave solutions 
and solid lines give results of the potential model calculations using AVig. The relative np energy 
is 1 MeV. 




(a) 3 Po admixture in 1 S'o (b) 1 Sq admixture in 3 Pq 

Figure 4: PV admixtures for the scattering states in Fig. 6L obtained by solving the inhomogeneous 
differential equations. 
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(a) l P\ admixture (b) 3 Pi admixture 

Figure 5: PV admixtures of deuteron: dashed lines are results using the plane wave Green's function 
and solid lines are calculations using AVis- 

but slightly differ in magnitudes. For the deuteron mixing, the l P\ state is induced only 
by p and u exchanges (AT = 0), while the 3 Pi is induced dominantly by the n exchange; 
both results are plotted in Fig. j|. Also shown by the dotted lines in the same figure are the 
solutions of the plane wave Green's function. Though these curves are similar in shape, the 
potential model calculation gives smaller amplitudes and different small-r radial dependence 
than plane waves. From now on, we only present results from the potential model calculation 
which is more realistic. 

As the impact of G$ on A LR is more important at backward angles, we first ex- 
amine the extreme case: 6 = 180°. Subsequently, we present results relevant to SAMPLE 
kinematics. The maximum Q 2 we consider is 0.15 (GeV/c) 2 , and the saturation behavior 
shown in Fig. || justifies truncation of the sum over final scattering states at total angular 
momenta J/ < 7. 

Fig. |7| indicates how the backward angle cross section and asymmetry vary with Q 2 , 
ranging from 0.01 to 0.15 (GeV/c) 2 at the QE peak. It is clear that the asymmetry due 
to two-body hadronic PV (plotted with a magnification of 100) is insignificant compared 
with the contribution of tree-level Z°-exchange plus radiative corrections, which includes 
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(a) Cross section 



(b) Asymmetry 



Figure 6: The saturation behavior of total cross section and asymmetry as functions of J ma x, the 
maximum total angular momentum for final states being included. 




Q 2 (GevV) Q 2 (GevV) 

(a) Cross section (b) Asymmetry 



Figure 7: The total cross section and PV asymmetry versus Q 2 , where the kinematics are con- 
strained to satisfy q 2 = 2m N uj and 6 = 180°. Note in (b), the hadronic PV contribution is 
multiplied by 100; it does not actually cross the Z° term in this Q 2 range. 

nucleon anapole effects. The curve for Z°-exchange asymmetry, plotted using the static 
approximation result in Ref. Q with parametrized nucleon form factors, shows the expected 
proportionality to Q 2 . The curve for hadronic PV shows a 0.05% correction to A^J at 
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Q 2 = 0.1 and a 0.3% correction at Q 2 = 0.04. We note that these results are consistent with 
simple scaling arguments as Eq. (|30"D. Although there is some enhancement for A^Jr as Q 2 
decreases, even at Q 2 ~ 0.01, near the threshold for QE kinematics, the correction is less 
than 5%. 

A detailed breakdown of various hadronic PV contributions is shown in Fig. ||. The 
deuteron mixing, rather insensitive to the Q 2 of the explored region, is the dominant con- 
tribution for Q 2 > 0.03. Its correction to AfJ at Q 2 = 0.1 is approximately 0.1%, and 
0.3% at Q 2 = 0.04. These values are consistent in the order of magnitude with the results 
of Ref. |l(J, which used Bonn potential to calculate the PC wave functions and the plane 
wave Green's function to calculate the parity mixture in deuteron. On the other hand, the 
final state mixing and PV meson exchange currents, though comparatively smaller, do have 
a combined contribution which could be as large as half of the contribution from deuteron 
mixing for Q 2 > 0.03. They are also more sensitive to Q 2 and become important when 
approaching the QE threshold. 

Away from the QE peak, the dependence of the cross section and asymmetry on final 
electron energy are shown in Fig. |9| for 194 and 120 MeV beams. Since the scattered 
electrons are detected via the Cerenkov radiation (the threshold is about 20 MeV), scattered 
electrons with E' e < 150 and 100 MeV, respectively, are detected. However, judging from 
the cross section plot, only regions about peak energy ± 20 and 10 MeV, respectively, are 
important for these two cases. When these asymmetries are further plotted as ratios to 
A LR , as shown in Fig. |10J, we observe that the correction could become as large as a few 
percent. Notice, however, that the corrections change sign roughly when crossing the QE 
ridge. Hence, corrections from these two regions cancel after integration, thereby, keeping 
the total correction small. A similar feature was also found in the calculation of Ref. fl24|, 
where PC two-body effects were considered. 

The setup of SAMPLE experiments actually cover the angular range from 130° to 170°, 



the average angles of the detectors are: 138.4°, 145.9°, 154.0°, and 160.5° |g5j. The corre- 



sponding cross sections and asymmetries are plotted in Fig. fLT| The general trend is that 
when the angle gets smaller, the cross section becomes larger and the asymmetry becomes 
smaller. However, the overall behaviors are not too different from the 9 = 180° case. 

Summarizing these observations, we conclude that the two-body hadronic PV effects in 
QE e-d scattering are negligible. However, the situation changes in the kinematic region of 
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Figure 8: The breakdown of various hadronic PV contributions to the asymmetry at QE kinematics, 
where D, F, and MEC refer to contributions from deuteron mixing, final state mixing, and PV meson 
exchange currents, respectively, and the solid line gives the total. 



threshold disintegration as shown in Fig. |12[ At Q 2 ~ lCr 4 (GeV/c) 2 , these two are com- 
parable 5 , and hadronic PV dominates when moving toward lower Q 2 region. Here again, 
the magnitude of Q 2 at which the hadronic PV and Z°-exchange contributions are com- 
mensurate is roughly what one would expect based on the simple scaling arguments of Eq. 

. The detailed breakdown given in Fig. shows that the final state mixing has the 
most important contribution and that PV meson exchange currents are also significant. The 
deuteron mixing, still rather independent of Q 2 evolution, becomes negligible. We also point 
out that while our calculation in this kinematic region is consistent with Hwang et a/.'s || 
at Q 2 > 0.0001 GeV 2 /c 2 , we obtain larger asymmetries as one approaches the threshold 
region. The reason is that we use a potential (AVis) which has a much softer core than the 
Reid soft-core potential. Thus, the behavior of the wave function at low energy and small 



5 The Z° asymmetry is plotted using the same formula from the static approximation. Although there are 
also two-body effects, the Q 2 -dependence still governs the overall behavior. 
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(a) Cross section 
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Figure 9: The total cross section and two-body hadronic PV asymmetry versus electron final energy, 
for 194 and 120 MeV incident beams and 180° scattering angle. The asterisk denotes the position 
of the QE peak. 




E c ' (MeV) E c ' (MeV) 

(a) 194 MeV beam (b) 120 MeV beam 

Figure 10: The ratio of asymmetry due to twodaody hadronic PV and Z°- exchange, versus final 
electron energy. Both kinematics are the same as in Fig. |j| 

distance is important for studies of hadronic PV at threshold, including experiments like the 
photo-disintegration of deuteron, radiative neutron capture, and neutron spin rotation. 
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E c ' (MeV) E.' (MeV) 



Figure 11: The same plots as Fig. || with four average detector angles of SAMPLE experiments: 
(a) 160.5°, (b) 154.0°, (c) 145.9°, and (d) 138.4°. The left panels are for the E e = 194 MeV and 
right panels for the E e = 120 MeV case. 

V. CONCLUSIONS 

The theoretical analysis of PV electron scattering asymmetries requires that one take into 
account effects which may, in principle, cloud the intended interpretation of an experimental 
result. In this study, we have analyzed the effects of parity violating NN interactions which 
give rise to a non- vanishing inelastic eD asymmetry at the photon point. Our results indicate 
that for the QE kinematics relevant to the SAMPLE experiment, these effects generate a 
negligible contribution to the PV asymmetry. Moreover, contributions arising from each side 
of the QE peak produce cancellations when integrated over detector acceptances, thereby 
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Figure 12: The asymmetries due to hadronic PV and Z°exchange versus Q 2 in the threshold 
electro-disintegration region, with small, fixed np relative energies. 

generating an additional suppression of the nuclear PV contamination. From this stand- 
point, then, the PV QE asymmetry provides a theoretically clean environment for studying 
electroweak nucleon form factors, such as G ( f > ^ T ~ l \Q 2 ). On the other hand, PV effects in 
the threshold region can become dominant, with asymmetries as large as a few xO.l ppm. 
Hence, near-threshold electro- or photo-disintegration of the deuteron could provide a tool 
for probing the PV NN interaction. 
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